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Integrals of the form 



Q,(&)« pll-fc 2 cos 4>]- ( i+»d4> 



where <&<1, andj is a positive integer, occur in a radiation field problem. Expressions 
for Q (k) and tti(k) have been derived in terms of complete elliptic integrals of the first 
and second kinds. Using these values, and the recursion formula 

(2j- 1) a-k^i(k) =40- l)0/-i(*) - (2j-3)Q,- 2 « 

12,-(A0 can be found for all values of j and k. A number of useful series expansions and other 
relations are given for 0,(ft), and tables are included for 0<j(l)<9 and 0<& 2 (0.01) <0.99. 

1. Introduction 

In calculating the radiation field off-axis from a uniform circular disk radiating according 
to an arbitrary angular distribution law, it has been convenient [l] 3 to use the method of 
Legendre polynomial expansion of the source-detector geometry [2]. Evaluation of the 
coefficients in these expansions leads to a family of integrals 

Q t <Jc)m p[l-F cos 4>]' u ^d4 (1) 

where 0<&<1. This function has a singularity at k^l which corresponds, in the problem 
described above, to the detector at the edge of the disk. 

The integral (1) is closely related to the complete elliptic integrals of the first and second 
kind. Thus it can be shown readily that [3, eq 291.00] 

Go(fc) = (VM)#M (2) 

where 

K a =2i 2 /(l+fc 2 ) (3) 

and K(k) is the complete elliptic integral of the first kind [3, eq 110.06] 

K(k)= r /2 (l-/c 2 sin 2 (A)^^. (4) 

Similarly, Qi(k) can be expressed [3, eq 291.01] as 

1 (fc)=(V2^)[l/(l~P)]J5;W (5) 

where [3, eq 110.07] 

E{k) m r 2 (1 -k 2 sin 2 0)+*efy. (6) 
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Uj(k) for j> 2 can also be developed in terms of these functions. Thus, using the inte- 
gration by parts formula given by Peirce-Foster [4, eq 317], it follows at once that 

(2j-l)(l-k")^k)=Mj~l)^-i(k)-(2j-3)n j - 2 (k). (7) 

This recursion formula is exact, and by using it and values of K(k) and E(k), toj(k) 
can be obtained for all values of j. For example, it is found by inserting (2) and (5) in (7) that 

^ 2 (/:)-(V2,A)[(l/3)(l-F)][4£W/(l-F)-^W] (8) 

and 

^(&)=(V2k/*)[(i/^ (9) 

It is clear that continuing in the same way 2j(k) can be obtained for all values of j in terms of 
K and E. Extensive tables of these complete elliptic integrals exist [3, 5, and 6] in addition 
to excellent polynomial approximations suitable for use on a high speed computer, due to 
Hastings [7], and a large number of exact and approximate analytical formulations [6, 8]. 
For values of k 2 approaching unity, the expansions (see, e.g., [5, eqs 773.3 and 774.3]) 

I£00=lri (4//) + (l/4){ln (4//)-1}k' 2 + . . . (10) 

and 

SW = l + (l/2){ln(4//)-l/2}/ 2 + . . ., (11) 



where k =^l — k 2 , can be used to compute ti (k) and &i(k) to high precision. 

This technique for finding ttj(k) has several serious disadvantages, however; e.g., (a) to 
find fi y from (7), particularly for large values of j, the precision of the Qj thus computed may 
be poor due to the buildup of round-off errors ; (b) in order to find Qj for a given k, it is first 
required that i2 , £l ly ti 2 , ■ • • &/~i be known. If only a few values of £lj(k) are required, this 
method may be painfully long and complicated. For these and similar reasons, alternate 
methods for calculating Q.j(k) are considered below. 

As will also appear in the sections which follow, the quantities tij(k) are related to an 
interesting collection of transcendental functions. In addition to the elliptic integrals 
above, it will be shown that Qj(k) is expressible in terms of Bessel functions with a purely 
imaginary argument, and in terms of error functions. The analysis of these functions thus 
involves a fascinating excursion into a number of curious bypaths in higher mathematical 
analysis, in addition to the practical application to the physical problem mentioned in the 
opening sentences above. 

2. Series and Other Expansions for Qy(fc) 

2.1. Simple Power Series j>0 

For small values of k, £lj(k) can be derived by a binomial theorem expansion of 
[1 — Pcos <t>]~ u+ V in powers of {¥ cos <j>), i.e., 

°>(*)=g (T*) C-^'J^COB^rf* (1) 



where 



Q^nl/i\(n-i)\ 



(2) 



is the familiar notation for the binomial coefficient 4 [9, par. 1.51], [10]. The integral in (1) can 
be shown to be [4, eq 498] 



t As indicated in [9] and [10], binomial coefficients of the form (J 1 ^ ) are easily generated using the recursions: 



£ cos<*<fo=[l + (-l)%/?/2) \\ 2 ( (3) 

and so vanishes when i is odd. Using this relation and the well known properties of the bi- 
nomial coefficients, it follows that 

Qj(k) = ±a m (j)k im (4) 

TO=0 

where 

a ^J) ( 64 )m (2^)! (2m+j)! VW W 

Equation (4) may easily be shown to be convergent for all values of j and lor &< 1. For 
k close to unity, however, it tends to converge rather slowiy; and for this reason, this power 
series formulation is most useful for k<C<^l. 

2.2. Difference Sum, ;/>() 

The behavior of (4) may be improved by using a technique previously developed for 
improving the convergence of series [8]. By this method, the series (4) is written in the iden- 
tically equivalent form 

00 CO 00 

n/*)=S«»(i)* 4 "=S«i(i)* 4 "+ Z>«( J)-<U)\ic A '" (6) 

TO=0 TO=0 772=0 

where a* (j) is defined, first, by the condition 

lim[a*(i)K(i)]=l (7) 

77l-»a> 

that is, a* is asymptotically equivalent to a m . From the infinite variety of sets of a* values 
satisfying this condition, a set is chosen such that the first sum on the right of eq (6) can be 
evaluated in closed form. If this technique is properly applied, the residual sum of the differ- 
ences [oi m (j) — a* n (j)W m converges much more rapidly than the original sum in a m (j)¥ m . 

For the present problem, it is shown in appendix A that by substituting Stirling's formula 
[11, p. 5] for the factorials in eq (5) such a set of a*(j)'s can be generated, of the form 



(8) 



where 



and 



ao)-^{§ 5 „("' + t i - 1 )+ s ii} 

,4,= (l/V2)2-"( 2 £) (9) 

B i} = (-l)-2-* @/[( 2 r f ) (£|.) ii-i)} (10) 

The notation ^ means that the finite sum over i vanishes when j=0. 

i>0 

Using (8), sums over a^(j)k im can then be identified with the closed-form expressions 

m=0 V. !>0 } 

From (8) we have, for j=0, 1 and 2 

a*(0) = (l/V2) ^, (12) 
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a* n (1) = (1/V2) 1 4+ (1/4) ^L. | , (13) 

and 

«*(2) = (1/V2){(16/3)(to+1)-1 + (3/32) ^jV (14) 

The corresponding sums from (11) are 

± a*(0)^=(l/V2)(-l#*) In (1-k*), (15) 

m=0 

± «*(l)^=(l/ A / 2){4/(l-^) + (l/4)(-l/^) In (1-f)}, (16) 

7W = 

and 

S a*(2)^=(l/ A ,/2){(16/3)/(l-F) 2 -l/(l-^) + (3/32)(-l/^) l n (1-& 4 )}. (17) 

7tt = 

Using (5) and the particular form of a* (j) given in (8), the residual sum in (6), viz, 

S=2 [ab(j)-o*(i)]iH- 
converges as fast as 23 [l/(ra + l) 2 ]& 4 ™. This may be contrasted with the simple power series 

771=0 

(4) which only converges as fast as V} m j ~ l k 4m . 

771=0 

An additional bonus is noted for the higher values of j. While X R for any given j varies 
by less than a factor of two over the entire range 0<k 2 < 1, the total value of the function Qj(k) 
varies from -k up to infinity over this same range. Thus, in the region of weakest convergence 
of 2#, i.e., 0.90 <P<1, the contribution of X R to fy(&) may be neglected in many cases. 

Also, in the right hand side of (11), the successive terms in the finite sum decrease very 
rapidly with increasing i, and the log-term, (— l/k 4 ) ln (1— k 4 ), is, for jf>l, 0<& 2 <1, of the 
same order of magnitude or less than the final term (i=j—l) in the finite sum. Thus, for 
example, for ,7>6, & 2 >0.90, values of Qj(k) accurate to one part in 10 8 may be computed by 
omitting from (11) the finite series terms for i>o, the log-term (—l/k 4 ) ln (1 —k 4 ), and the entire 
sum 2 R on the right side of (6). 

This ease of attaining toj(k) values of high precision for larger' proved useful in evaluating 
disk-source Legendre coefficients [1] for large /, which involve smaller and smaller differences 
between larger and larger numbers. 

2.3. tij(k) in Terms of Bessel Functions, j>0 

The differencing technique illustrated for the sum in section 2.2 above may often profitably 
be applied to integrals as well, i.e., by writing 

f/(*)<fo= [\f(st>)d4>+ fV(*) -/*(«)]<&• (18) 

Ja Ja Ja 

In this expression /*(<£) is a function, integrable in closed form over the region a <<t><b, which 
is a suitable approximation to /(<£). Applied to ttj(k), it will be noted from (1 — 1) that 

/(*) = (l-Ar 2 cos*)-» (19) 

where 

n=j+h (20) 
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and 

a = 0, b = ir. (21) 

For this case, it is convenient to take 

f*(cj)) = e +n(k2coa + ) + (;n/2)(k 2 cos 0) V^ 27 ^* 2 ' 08 ^, (22) 

an expression derived by retaining terms through (k 2 cos <£) 3 in the expansion of j(4>) as given 
by (19). Then 

f V*(*)<fy= f'exp [+ra 2 (Jfc 2 cos 0)]rf«+ (n/2)k 2 f * (cos 2 0) exp |~(W|)(& 2 cos </>)!#. (23) 

Xow from the expression [12, p. 79, Xo. 9] 

/>(g)s r( y +^rft) Jo^ (sin2 0)y[d:g cos *!''*> (24) 

where 7,(2) is the Bessel function of order ^ with a purely imaginary argument, i.e., 

L{z)=i' v J v {iz). (25) 

From (24) 

I (z) = (1Itt) I "exp [±zcos<t>](kt> (26a) 

Jo 

and 

Ti(z) = (z/t) I * (sin 2 0) exp [±z cos 0]d0. (26b) 

Jo 

Using (26) to evaluate the integrals in (23), it is found that 
(l/7r)G#) = +J [F(^ 

+ (l/^)| , o T [(l-^ 2 cos0)-^^-exp{ + (i+i)(Fcos0)} 

- (!) (i+i) (£ 2 cos 0) 2 exp { O'+t) (^ cos *)}+.. . ]*/>. (27) 

By expanding 

(l—x)- n —e +nx —inx 2 e+( n + l)x 

in a power series in x=k 2 cos <j> and integrating term by term, the integral in (27) may be shown 
to be equal to 

S(l/m!) 2 [+ ( ^m)\ ' (^]— 4 m (2i+l) 2m -18(2/3) 2m • m(2m-l)(2j+l)(6j+7) 2 -- 2 ](*764)»' 

= (3/64) (i+|) (*+§)*»+.... (28) 

Further manipulation and simplification in the case of (28) is possible, but these procedures 
have not proven particularly rewarding. The greatest importance of eq (27), which it will 
be observed is exact, arises from the relation shown by this equation between Qj(k) and the 
Bessel functions. For numerical calculation, forms other than (27) and (28) may be more 
convenient, except possibly for some special cases. 

It will be noted that /*(<£) as given by (22) is a good approximation to /(</>) only for small 
values of k 2 , since it is derived from the first three terms of the power series expansion of (19). 
Consequently, (27) is essentially useful only for |&|<<Q. To realize the full potentialities of 
the formulation (18), it would be necessary to find a form for /*(<£) valid over the whole range 
0<C&<T1 and integrable in closed form. Such an approximation is explored in section 2.5 
below. 



2.4. ttj(k) in Terms of Error Functions (the Method of Steepest Descent) 
By writing 

Jo 
so that 

*■(*) = (/+*) to (l-Fcos</>) (29) 

it is possible to obtain another interesting form for the function. Expanding F(cj>) in a 
Maclaurin series [9, par. 6.104], it can readily be shown that all the odd derivatives of F(<f>) 
vanish at the origin and that 

^(k)=e- Fm f%- (V ' ) V[i-(i/4!)jW)(o)^--(l/6!)F (W (0)^ l + . . .]d</> 

Jo 
where 

(X/7r) 2 =(l/2)F (2) (0). (30) 

Inserting the value of 2^(0) and making the transformation t=(\/7r)<t), it follows at once that 
a#) = (7r/X)(l-P)-W>[^ . . . ] (31) 

where 

L p (x)= (" t*'e-'*dt. (32) 

Jo 

The integrals L P (x) can be expressed in terms of the incomplete gamma [13] or factorial [11] 
functions, L p (x)=%- (p— J,# 2 )! but an alternate procedure is perhaps more informative. Inte- 
grating (32) by parts, from the resulting recursion formula 

L p+1 (x) = (l/2)\(2p+l)L p (x)-x 2 "+ l e-^] (33) 

the values of L p (x) can all be found in terms of L (x). But this is simply 

Lo(z) = W*l2)evtz (34) 

where the error function [5, par. 5.90] is 

erfz = (2/V^) p e- %2 dt. 
By this method then, it is found that 

!*(*) = (3/4) [L (x)-xe~* 2 (l+| x 2 )] (35a) 

Zb(s) = (15/8) [Lo(z)-ze~* (l+| x 2 +^ x*jj- (35b) 

The higher derivatives of F(<t>) at the origin, F (p) (0), are readily determined by expanding 
F(4>) in a power series and recognizing that F (p) (0) is p! times the coefficient of <j> p . In this 
way, it is found that 

F^(0) = + (j+±)k*(l-k>)- 1 (36a) 

^ (4) (0) = -(i+4)^ 2 (l-^ 2 )" 2 [2/: 2 + l] (36b) 

2T<«)(0) = + y+J)jfc»(l-jfc«)- 8 [l6fc 4 +13t I +l]. (36c) 

From (36 A) and (30), 



X=(7r/2)V'(2i+l)P/(l-P) (37) 
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and combining this and the previous results with (31), it follows that 

Q#WW(2j+l)F(l-F)-'[+erfX+(l/2)(2^ 

-(l/12)(2i+l)- 2 (io+P+p){erf \-(2/^)(\e^ 2 ) (l+§ X 2 +~ \')\ + . . .] . (38) 

In using these relations for small k 2 or X values, the problem of small differences between large 
quantities is resolved by the series 

{erf X-(2/V^)(\6- x2 )(l + |X 2 )}-(2/V^) (4X 5 /15) {1- (5/7)X 2 +(5/18)X^-(5/66)X 6 

+ (5/312) X 8 - (1/360) X 10 + (1/2448) X 12 + . . .} (39) 
and 

^erf X-(2/V?)(Xe-* 2 ) (l+| X 2 +^| X 4 )\ 

= (2/V^)(8X 7 /105) {l + (77/16)X 2 +(7/22)X 4 -(7/78)X 6 +(7/360)X 8 -(7/2040)X 10 + . . . } . (40) 

For large X values, i.e., large j's, erf X and the bracket terms involving this quantity are all 
close to unity, and 

lim Q J W = VW(2i+DF(l-F)-^[l + (l/2)(2i+l)- 1 {l+^| 

-MflrrD-{l+ 3 W ! }+...]. (41) 

2.5. Numerical Integration Techniques for ilj(k) 
To apply eq (18) in a more general way than was used in section 2.3 above, let 

/*(<£) = (l-F)-« (42) 

where n=j+Y 2 , as before (eq 20). Then from (18) 



Qj(k) = 7r(l-k 2 )- n + (* [(1-k 2 COS<t>)- n -(l-k 2 )-"]< 



/*. 



It will be observed that the integrand here contains no singularities, even at the origin for 
jfc=l. With the change of variable <t>=irx } this becomes 

9.j{k)=7r(l-k 2 )- n +7r C [{1-k 2 cos (wx)- fl }-(l-k 2 )- n }dx. 

This integral is now in standard form for the use of numerical integration techniques. By 
fitting a polynomial of the JV'th degree to E(^), it can be shown that 

f 5(x)dx^f] w N (i)E(i/N) If] w N (i) (43) 

,7 i=0 / i=0 

where the w N (i)'s are constants (given for example in the introduction to the Tables of La- 
grangian Interpolation Coefficients [14] where the integrals are expressed in slightly different 
form). It should be noted that (43) can never be expected to give exact answers; but by in- 
creasing the value of N, closer and closer approximations to ; (&) can be obtained over the 
whole range of k, although this procedure does not necessarily result in improved accuracy at 
any given value of k. In the present case, where for certain values of N, cos (iw/N) can be 
expressed in simple closed form, this technique leads to an interesting series of approximations 



for ttj(k) as given below. Note that all these results yield exactly the right answers at 
k=0, Qj(0) = Tr 

.,1 



N=2 
Q#) = (+6)[+4+(l-F)-»+(l+& 2 )-"] (44) 



iV=3 



Bj(i) = (ir/8) [+(l-F)-*+(l+P)-"+3(l-iF)~ K +3(l+^p)~" (45) 

iV=4 
ftXt) = (x/90) r+12+7(l-F)- B +7(l+F)-"+32('l-4Fy"+32 A+-Ljfe»Y" (46) 

iV=5 
B i (Jfc)=(x/288) r+19(l-Ar 2 )-"+19(l+t 2 )-' ! +75 A-^tl jfciV" 

+75 (l+fe f )""+50 (l-^T 1 ^)""+60 (l+^ * 2 )~" (47) 

N=6 
fi#) = (W840)[ + 272+41(l-P)-"+41(l+P)-"+216(l-^/3P)-' ! 

+216(1+W3& 2 )~ n +27(1— §& 2 )-"+27(l+ffc 2 )- n ] (48) 

7V=10 
Q,(jfc) = (ir/598,752) +427,368+16,067(1 -& 2 )- n +16,067(l+& 2 )-" 

+ 106,300 h-yfc^jpj +106,300 (i+^3*p\ -48,525 (l- 2 ^ *•)"" 

-48,525(1+^1^) "+272,400 (i-^t^fapj +272,400 (l+y ^=^ A» J 

-260,550 (l— fc^ifc 2 ) -260,550 N+^ifr 2 ) • (49) 



The interesting and beautiful symmetries in these approximations have been emphasized in 
the formulas above. It is evident that by using a more complicated expression for /*(</>) than 
(42), better agreement can be achieved. Also, by going to larger values of N, or by dividing 
the range of <j> into smaller intervals, formulas of greater accuracy can be developed. For large 
values of j and k 2 , or their product 

lim Q i (k)=w(l-k i )- (i+ i'w N (0) /£; w N (i) 

#!->■= / 1=0 
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or inserting the expression from w N (i) given by Whittaker and Robinson [15] 



lim fi ; (*)=x(l-A 2 )- (y+ * 



N 



Jo \ N 



dt 



(50) 



where, it will turn out, since the w N (i)'s are all integers, the last term here is a rational fraction. 

3. Numerical 



It is immediately apparent from the definition of the function that 

Qj(0) = t and 0,(1 )=«. 



(1) 



The intermediate behavior of ttj(k) versus k 2 is displayed for <.7*(1) <9 as a family of curves 
in figure 1. These curves are based on values of fy(&) listed for 0<j(l) <9 and <& 2 (0.01) < 
0.99 in table 1. 

Table 1 was computed in "double precision" on the National Bureau of Standards IBM 
7090 electronic computing* machine, using the difference-sum (2.2-6) with the particular 
a m(j) given in (2.2-8). The results are given t o eight digits, which were rounded off from a table 
correct to at least nine digits throughout, so all eight digits can be considered significant. 
The results in the region of poorest convergence (high k, low j) of the difference sum were 
verified by ;i ten-digit hand calculation using the recursion expression (1 7) and the high-Jfc 
elliptic integral expansions (1—10) and (1-11). 

The lengthy developments in section 2 can be considered to yield various approximations 
to the value of i}j(k). In tables 2 to 6 numerical examples are presented which will give some 
idea of the usefulness of the techniques. For each example, the values ^"=0 and 8, and k 2 = 
0.01 and 0.99 are used for purposes of illustration. 




Figure 1. Family o/12,(k) curves based on numerical 
results in table 1. 



667805— G3 2 



Table 1.* Values of fij(k) for 0<j(l) <9 onrf 0<k 2 (0.01)<1.00 

These were computed usine the difference sum (2.2—6) with the a m *(j) in (2.2—8). A tolerance of 
<10-u was set on the fractional contribution of the first neglected series term [a m (f)—a m *(j)]k* t * 
which resulted in nine or more significant figures throughout, here rounded off to eight 
significant figures. 
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3. 1532035 


[0) 


3. 2004553 


(0) 
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3. 3950431 


(0) 


3.5463805 (0) 


.15 


3. 1550119 


,0) 


3. 2093019 
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3. 3017684 
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3. 4344882 
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3.6105249 (0) 
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3. 1568872 


:o) 
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3. 4772340 
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r 0) 


3. 2518347 


(0) 


3. 4042430 
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(0) 


3. 4344610 


(0) 


3. 6839559 


(0) 


4.0230902 (0) 


.21 


3. 1682158 


0) 


3. 2773746 


(0) 


3. 4666396 


(0) 


3. 7453695 


(0) 


4.1264381 (0) 


.22 


3. 1708833 


0) 


3. 2912770 


(0) 


3. 5008493 


(0) 


3.8110474 


(0) 


4.2377233 (0) 


.23 


3. 1736890 


,0) 


3. 3059591 


(0) 


3.5371053 


(0) 


3. 8812046 


(0) 


4.3574565 (0) 


.24 


3. 1766349 


0) 


3.3214413 


(0) 


3. 5756686 


(0) 


3. 9560749 


(0) 


4.4861987 (0) 


.25 


3. 1797233 


,0) 


3. 3377450 


(0) 


3. 0104408 


(0) 


4.0359125 


(0) 


4.6245662 (0) 


.26 


3. 1829566 


0) 


3. 3548935 


(0) 


3. 6595943 


(0) 


4. 1209941 


(0) 


4.7732360 (0) 


.27 


3. 1863371 


0) 


3.3729113 


(0) 


3. 7052127 


(0) 


4.2116207 


(0) 


4.9329530 (0) 


.28 


3. 1898076 


0) 


3. 3918249 


(0) 


3.7534115 


(0) 


4.3081200 


(0) 


5. 1045300 (0) 


.29 


3. 1935508 


0) 


3. 4116624 


(0) 


3. 8043080 


(0) 


4. 4108487 


(0) 


5.2888871 (0) 


.30 


3. 1973897 


0) 


3. 4324538 


(0) 


3. 8580313 


(0) 


4. 5201954 


(0) 


5.4869998 (0) 


.31 


3. 2013874 


0) 


3. 4542312 


(0) 


3. 9147107 


(0) 


4. 6365838 


(0) 


5.6999702 (0) 


.32 


3. 2055471 


0) 


3. 4770287 


(0) 


3. 9745127 


(0) 


4. 7004758 


(0) 


5.9290085 (0) 


.33 


3. 2098724 


0) 


3. 5008829 


(0) 


4. 0375787 


(0) 


4. 8923759 


(0) 


6.1754521 (0) 


.34 


3. 2143669 


0) 


3. 525S327 


(0) 


4. 1040871 


(0) 


5. 032N355 


(0) 


6.4407815 (0) 


.35 


3. 2190345 


0) 


3. 5519197 


(0) 


4. 1742238 


(0) 


5. 1824573 


(0) 


6.7266372 (0) 


.36 


3. 2238793 


0) 


3. 5791883 


(0) 


4. 2481901 


(0) 


5. 3419010 


(0) 


7.0348393 (0) 


.37 


3. 2289056 


,0) 


3. 0076802 


(0) 


4.3202030. 


(0) 


5.5118921 


(0) 


7.3074111 (0) 


.38 


3.2341182 


0) 


3. 0374640 


(0) 


4.4085000 


(0) 


5.0932230 


(0) 


7.7206040 (0) 


.39 


3. 2395217 


,0) 


3. 0085702 


(0) 


4. 4953348 


(0) 


5. 8807072 


(0) 


8. 1149281 (0) 


.40 


3. 2451213 


0) 


3. 7010809 


(0) 


4. 5809850 


(0) 


6. 0934851 


(0) 


8.5351804 (0) 


.41 


3. 2509225 


0) 


3. 7350405 


(0) 


4. 6837518 


(0) 


6. 3144352 


(0) 


8.9905145 (0) 


.42 


3. 2569310 


0) 


3. 7705219 


(0) 


4. 7859020. 


(0) 


6. 5507857 


(0) 


9.4844263 (0) 


.43 


3.2631530 


0) 


3. 8075968 


(0) 


4. 893973S 


(0) 


6. 8038277 


(0) 


1.0020808 (1) 


.44 


3. 2095947 


0) 


3. 840)3422 


(0) 


5. 0081739 


(0) 


7.0749911 


(0) 


1.0604277 (1) 


.45 


3. 2762632 


,0) 


3. 8868412 


(0) 


5. 1289871 


(0) 


7. 3658616 


(0) 


1. 1239660 (1) 


.46 


3. 2831656 


0) 


3. 9291828 


(0) 


5. 2508770 


(0) 


7. 6782009 


(0) 


1. 1932671 (1) 


.47 


3. 2903097 


(0) 


3. 9734032 


(0) 


5. 3923515 


(0) 


8. 0139705 


(0) 


1.2689713 (1) 


.48 


3. 2977036 


[0) 


4. 0197800 


(0) 


5. 5359670 


(0) 


8. 3753580 


(0) 


1.3518054 (1) 


.49 


3. 3053561 


'0) 


4. 0082033 


(0) 


5. 6883349 


(0) 


8. 7048084 


(0) 


1.4425903 (1) 


.50 


3. 3132763 


(0) 


4. 1190101 


(0) 


5. 8501280 


(0) 


9. 1850001 


(0) 


1.5422873 (1) 


.51 


3. 3214743 


:o) 


4. 1721755 


(0) 


0. 0220875 


(0) 


9. 6391874 


(0) 


1.6519570 (1) 


.52 


3. 3299604 


(0) 


4. 2278838 


(0) 


6.2050323 


(0) 


1. 0130049 


(1) 


1.7728428 (1) 


.53 


3. 3387459 


[0) 


4. 2802954 


(0) 


0. 3998082 


(0) 


1. 0663340 


(1) 


1.9003672 (1) 


.54 


3. 3478428 


[0) 


4. 3475784 


(0) 


6. 6075997 


(0) 


1. 1241689 


(1) 


2.0541717 (1) 


.55 


3. 3572639 


(0) 


4.4119159 


(0) 


6. 8293426 


(0) 


1. 1870679 


(1) 


2.2181554 (1) 


.56 


3. 3670230 


(0) 


4. 4795070 


(0) 


7. 0063401 


(0) 


1. 2556002 


(1) 


2.4005228 (1) 


.57 


3. 3771348 


:o) 


4. 5505724 


(0) 


7. 3199790 


(0) 


1.3304138 


(1) 


2.6038421 (1) 


.58 


3. 3876151 


:o) 


4. 0253495 


(0) 


7. 5918137 


(0) 


1.4122502 


(1) 


2.8311150 (1) 


.59 


3. 3984809 


(0) 


4. 7041020 


(0) 


7. 8835849 


(0) 


1. 5019596 


(1) 


3.0858639 (1) 


.60 


3. 4097506 


[0) 


4. 7871193 


(0) 


8. 1972534 


(0) 


1. 0005209 


(1) 


3.3722376 (1) 


.61 


3. 4214439 


(0) 


4. 8747204 


(0) 


8. 5350309 


(0) 


1. 7090647 


(1) 


3.6951427 (1) 


.62 


3. 4336822 


(0) 


4. 9072582 


(0) 


8.8994209 


(0) 


1. 8289010) 


(1) 


4.0604078 (1) 


.63 


3. 4461887 


(0) 


5. 0651235 


(0) 


9. 2932050 


(0) 


1. 9015571 


(1) 


4.4749869 (1) 


.64 


3. 4592883 


(0) 


5. 1087509 


(0) 


9.7198032 


(0) 


2. 1088134 


(1) 


4.9472173 (1) 


.65 


3. 4729085 


[0) 


5. 2786244 


(0) 


1. 0182736 


(1) 


2. 2727023 


(1) 


5.4871440 (1) 


.66 


3. 4870789 


(0) 


5. 3952852 


(0) 


1.0686311 


(1) 


2. 4558693 


(1) 


6.1069344 (1) 


.67 


3. 5018320 


^0) 


5. 5193397 


(0) 


1. 1235424 


(1) 


2. 6610551 


(1) 


0.8214090 (1) 


.68 


3. 5172034 


[0) 


5. 6514705 


(0) 


1. 1835738 


(1) 


2. 8917967 


(1) 


7.0487274 (1) 


.69 


3. 5332319 


:o) 


5. 7924479 


(0) 


1. 2493S34 


(1) 


3. 1522555 


(1) 


8.6112803 (1) 


.70 


3. 5499606 


0) 


5. 9431449 


(0) 


1. 3217398 


(1) 


3.4474411 


(1) 


9.7308612 (1) 



See footnote at end of table, p. 12. 
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Table 1.* Values o/fij(k) for 0<j(l) <9 and 0<k 2 (0.01)< 1.00— Continued 



A-2 


fio 


ftl 


tt 2 


fi 3 


n 4 


.71 


3.5674369 (0) 


6. 1046647 (0) 


1.4015448 (1) 


3.7834208 (1) 


1. 1060216 (2) 


.72 


3.5857131 (0) 


6.2778123 (()) 


1.4898627 (1) 


4. 1675905 (1) 


1.2625110 (2) 


.73 


3.6048476 (0) 


6. 4642206 CO) 


1.5879565 (1) 


4.6090284 (1) 


1.4487127 (2) 


.74 


3.6249055 (0) 


6. 6662823 (0) 


1.6973345 (1) 


5. 1189617 (1) 


1.6717465 (2) 


.75 


3.6459595 (0) 


6.8827406 (0) 


1.8198097 (1) 


5.7113855 (1) 


1.9408189 (2) 


.76 


3.6680915 (0) 


7. list 1284 (0) 


1.9575775 (1) 


6.4038975 (1) 


2.2679546 (2) 


.77 


3.6913941 (0) 


7.3753315 (0) 


2. 1133163 (1) 


7.2188313 (1) 


2.6690316 (2) 


.78 


3.7159722 (0) 


7.6556690 (0) 


2.2903221 (1) 


8. 1848192 (1) 


3.1652641 (2) 


.79 


3.7419457 (0) 


7.9(129956 (0) 


2.4926875 (1) 


9.3389738 (1) 


3.7853627 (2) 


.80 


3.7694524 (0) 


8.3013360 (0) 


2.7255455 (1) 


1.0729980 (2) 


4.5687313 (2) 


.81 


3.7986513 (0) 


8.6755606 (0) 


2.9954048 (1) 


1.2422547 (2) 


5.5702880 (2) 


.82 


3.8297278 (0) 


9.0916202 (0) 


3.3106179 (1) 


1.4503942 (2) 


6.8678797 (2) 


.83 


3.8628993 (0) 


9.5568633 (0) 


3.6820480 (1) 


1.7093748 (2) 


8.5739421 (2) 


.84 


3.8984232 (0) 


1.0080473 (1) 


4.1240339 (1) 


2.0358783 (2) 


1.0854291 (3) 


.85 


3.9366074 (0) 


1.0674075 (1) 


4.6558190 (1) 


2.4536454 (2) 


1.3959246 (3) 


.86 


3.9778240 (0) 


1. 1352615 (1) 


5.3037170 (1) 


2.9972313 (2) 


1.8276781 (3) 


.87 


4.0225290 (0) 


1.2135628 (1) 


6. 1044816 (1) 


3.7182365 (2) 


2.4426513 (3) 


.88 


4.0712886 (0) 


1.3049152 (1) 


7.1107150 (1) 


4.6960084 (2) 


3.3432588 (3) 


.89 


4. 1248185 (0) 


1.4128679 (1) 


8.3998552 (1) 


6.0567809 (2) 


4.7056582 (3) 


.90 


4. 1840408 (0) 


1.5423880 (1) 


1.0089734 (2) 


8.0095478 (2) 


6.8473464 (3) 


.91 


4.2501717 (0) 


1.7006477 (1) 


1.2366829 (2) 


1.0917125 (3) 


1.0373314 (4) 


.92 


4.3248613 (0) 


1.8984009 (1) 


1.5540620 (2) 


1.5446583 (3) 


1.6516804 (4) 


.93 


4.4104268 (0) 


2. 1525401 (1) 


2.0155731 (2) 


2.2914615 (3) 


2.8010734 (4) 


.94 


4.5102677 (0) 


2.4912083 (1) 


2.7244578 (2) 


3.6165464 (3) 


5.1590976 (4) 


.95 


4.6296635 (0) 


2.9650431 (1) 


3.8964807 (2) 


6.2117605 (3) 


1.0636320 (5) 


.96 


4.7774745 (0) 


3.6752779 (1) 


6.0473488 (2) 


1.2060267 (4) 


2.5819863 (5) 


.97 


4.9703403 (0) 


4.8580263 (1) 


1.01)79679 (3) 


2.8419636 (4) 


8. 1144742 (5) 


.98 


5.2456872 (0) 


7.2213223 (1) 


2.3872660 (3) 


9. 5301054 (4) 


4.0851238 (6) 


.99 


5.7231768 (0) 


1.4303674 (2) 


9.4878355 (3) 


7.5852838 (5) 


6.5002882 (7) 


1.00 


CO 


CO 


CO 


CO 


CO 



w 


&5 


ih 


n 7 


ih 


S2 9 


0.00 


7T 


IT 


IT 


■K 


7T 


.01 
.02 
.03 
.04 
.05 


3. 1444016 (0) 
3. 1528418 (0) 
3. 1669537 (0) 
3. 1868051 (0) 
3.2124917 (0) 


3. 1454234 (0) 
3, 1569389 (0) 
3. 1762090 (0) 
3.2033511 (0) 
3.2385310 (0) 


3.1466027 (0) 
3. 1616703 (0) 
3. 1869086 (0) 
3.2225075 (0) 
3.2687360 (0) 


3.1479395 (0) 
3. 1070377 (0) 
3. 199060S (0) 
3.2443006 (0) 
3.3031717 (0) 


3. 1494340 (0) 
3. 1730429 (0) 
3.2126747 (0) 
3.2687599 (0) 
3.3419109 (0) 


.06 
.07 
.08 
.09 
.10 


3.2441379 (0) 
3.2818977 (0) 
3.3259565 (0) 
3.3765321 (0) 
3.4338769 (0) 


3.2819649 (0) 
3.3339215 (0) 
3.39 17246 (0) 
3.4647569 (0) 
3.5444634 (0) 


3.3259458 (0) 
3.3945750 (0) 
3.4751540 (0) 
3.5683119 (0) 
3.6747845 (0) 


3.3762170 (0) 
3.4641155 (0) 
3.5676927 (0) 
3.6879319 (0) 
3.8259893 (0) 


3.4329324 (0) 
3.5428347 (0) 
3.6728507 (0) 
3.8244570 (0) 
3.9993993 (0) 


. 11 
. 12 
.13 

. 14 
.15 


3.4982801 (0) 
3.5700701 (0) 
3.0496177 (0) 
3. 7373392 (0) 
3.8337006 (0) 


3.6343569 (0) 
3.7350228 (0) 
3.8471259 (0) 
3.9714180 (0) 
4. 1087457 (0) 


3.7954235 (0) 
3.9312075 (0) 
4.0832553 (0) 
4.2528403 (0) 
4.4414078 (0) 


3.9832113 (0) 
4. 1611549 (0) 
4.3616121 (0) 
4.5866377 (0) 

4.8385832 (0) 


4. 1997232 (0) 
4. 427S102 (0) 
4.6864215 (0) 
1.9787483 (0) 

5.3084726 (0) 


.16 

!l8 

.19 
.20 


3.9392217 (0) 
4.0544809 (0) 
4. 1801214 (0) 
4.3168573 (0) 
4.4654806 (0) 


4.2600610 (0) 
4.4264317 (0) 
4.6090547 (0) 
4.8092705 (0) 
5.0285795 (0) 


4.6505957 (0) 
4.8822565 (0) 
5.1384852 (0) 
5.4216493 (0) 
5.7344248 (0) 


5. 1201345 (0) 
5.4343576 (0) 
5.7847509 (0) 
6. 1753062 (0) 
6.6105807 (0) 


5 6798373 (0) 
6.0977301 (0) 
6.5677814 (0) 
7.0964805 (0) 
7.6913118 (0) 


.21 
.22 
.23 
.24 
.25 


4.6268698 (0) 
4.8019991 (0) 
4.9919487 (0) 
5. 1979173 (0) 
5.4212351 (0) 


5.2686619 (0) 
5. 5313990 (0) 
5.8188984 (0) 
6. 1335234 (0) 
6.4779255 (0) 


6.0798376 (0) 
6.4613109 (0) 
6.8827211 (0) 
7.3484615 (0) 
7.8635180 (0) 


7.0957805 (0) 
7.6368586 (0) 
8.2406300 (0) 
8.9149065 (0) 
9.6686559 (0) 


8.3609171 (0) 
9.1152864 (0) 
9.9659856 (0) 
1.0926426 (1) 
1.2012186 (1) 


.26 
.27 
.28 
.29 
.30 


5.6633798 (0) 
5.9259941 (0) 
6.2109057 (0) 
6. 5201504 (0) 
6.8559986 (0) 


6.8550829 (0) 
7.2683447 (0) 
7.7214814 (0) 
8.2187441 (0) 
8.7649324 (0) 


8.4335554 (0) 
9.0650197 (0) 
9.7652570 (0) 
1.0542653 (1) 
1.1406796 (1) 


1.0512189 (1) 
1.1457381 (1) 
1.2517935 (1) 
1.3709692 (1) 
1.5051001 (1) 


1.3241396 (1) 
1.4635194 (1) 
1.6218282 (1) 
1.8019586 (1) 
2.0073050 (1) 


.31 
.32 
.33 
.34 
.35 


7.2209852 (0) 
7.6179448 (0) 
8. 0500508 (0) 
8.5208624 (0) 
9.0343771 (0) 


9.3654742 (0) 
1.0026517 (1) 
1.0755037 (1) 
1. 1558963 (1) 
1.2447325 (1) 


1.2368671 (1) 
1.3440887 (1) 
1.4637943 (1) 
1.5976553 (1) 
1.7476017 (1) 


1.6563161 (1) 
1.8270952 (1) 
2.0203264 (1) 
2.2393861 (1) 
2.4882297 (1) 


2.2418611 (1) 
2.5103359 (1) 
2.8182967 (1) 
3.1723419 (1) 
3.5803135 (1) 


.36 
.37 
.38 
.39 
.40 


9. 5950927 (0) 
1.0208078 (1) 
1.0879058 (1) 
1. 1614509 (1) 
1.2421771 (1) 


1.3430429 (1) 

1.4520056 (1) 
1.5729710 (1) 
1. 7074894 (1) 
1.8573456 (1) 


1.9158678 (1) 
2.1050462 (1) 
2.3181518 (1) 
2.5587002 (1) 
2.8308011 (1) 


2.7715029 (1) 
3.0946757 (1) 
3.4642066 (1) 
3.8877418 (1) 
4.3743600 (1) 


4.0515567 (1) 
4.5972404 (1) 
5.2307520 (1) 
5.9681868 (1) 
6.8289570 (1) 



See footnote at end of table, p. 12. 
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Table 1.* Values o/fij(k) for 0<j(l) <9 and 0<k 2 (0.01)< 1.00— Continued 



fe2 


fi 5 


fie 


S2 7 


fig 


fis 


.41 
.42 
.43 
.44 
.45 


1.3309183 (1) 
1.4286240 (1) 
1.5363771 (1) 
1.6554164 (1) 
1.7871617 (1) 


2.0245991 (1) 
2.2116313 (1) 
2.4212044 (1) 
2.6565293 (1) 
2.9213497 (1) 


3. 1392711 
3. 4897718 
3. 8889761 
4. 3447730 
4. 8665169 


(1) 
(1) 
(1) 
(1) 
(1) 


4. 9348723 
5. 5821923 
6.3317936 

7. 2022780 
8. 2160859 


(1) 
(1) 
(1) 
(1) 

(1) 


7.8365504 (1) 
9.0194823 (1) 
1.0412490 (2) 
1.2058044 (2) 
1.4008261 (2) 


.46 
.47 
.48 
.49 
.50 


1.9332450 (1) 
2.0955464 (1) 
2.2762382 (1) 
2.4778370 (1) 
2.7032673 (1) 


3.2200424 (1) 
3.5577398 (1) 
3.9404785 (1) 
4.3753812 (1) 
4.8708801 (1) 


5. 4653362 
6. 1545121 

6. 9499504 

7. 8707665 
8. 9400186 


(1) 
(1) 
(1) 
(1) 
(1) 


9 4003772 
1. 0788151 
1. 2419646 
1. 4344125 
1. 6622140 


(1) 
(2) 
(2) 
(2) 
(2) 


1.6327337 (2) 
1.9094667 (2) 
2.2408845 (2) 
2.6392830 (2) 
3.1200644 (2) 


.51 

.52 
.53 
.54 
.55 


2.9559372 (1) 
3.2398316 (1) 
3.5596245 (1) 
3.9208178 (1) 
4.3299101 (1) 


5.4369916 (1) 

6.0856567 (1) 
6.8311633 (1) 
7.6906742 (1) 
8.6848881 (1) 


1. 0185627 
1. 1641533 
1. 3349170 
1. 5359335 
1. 7734580 


(2) 
(2) 
(2) 
(2) 
(2) 


1. 9328437 
2. 2555083 
2.64192HS 
3. 100307.1 
3. 6670460 


(2) 
(2) 
(2) 
(2) 
(2) 


3.7026116 (2) 
4.4114329 (2) 
5.2776727 (2) 
6.3411131 (2) 

7.6528445 (2) 


.56 
.57 
.58 
.59 
.60 


4.7946069 (1) 
5.3240797 (1) 
5.9292896 (1) 
6.6233924 (1) 
7.4222490 (1) 


9.8388739 (1) 
1. 1183126 (2) 
1.2754912 (2) 
1.4599991 (2) 
1.6774836 (2) 


2. 0552302 
2. 3908727 

2. 79241(H) 
3. 2749532 

3. 8575879 


(2) 
(2) 
(2) 
(2) 
(2) 


4. 3469219 
5. 1751711 

6. 1890817 
7. 4365900 
8. 9797058 


(2) 
(2) 
(2) 
(2) 
(2) 


9.2788504 (2) 
1.1304851 (3) 
1.3842894 (3) 
1.7040380 (3) 
2.1092534 (3) 


.61 
.62 
.63 
.64 
.65 


8.3450710 (1) 
9.4152431 (1) 
1.0661378 (2) 
1.2118673 (2) 
1.3830682 (2) 


1.9349511 (2) 
2.2411423 (2) 
2.6070254 (2) 
3.0464479 (2) 
3.5770065 (2) 


4. 5645739 
5.4269316 

6. 4845852 

7. 7892704 
9. 4085123 


(2) 
(2) 
(2) 
(2) 

(2) 


1. 0899157 
1. 3300702 
1. 6324220 
2. 0155347 
2. 5043262 


(3) 
(3) 
(3) 
(3) 
(3) 


2.6259746 (3) 
3.2891909 (3) 
4.1462859 (3) 
5.2619563 (3) 
6.7253039 (3) 


.66 
.67 
.68 
.69 
.70 


1.5851615 (2) 
1.8249392 (2) 
2. 1109677 (2) 
2.4541303 (2) 
2.8683593 (2) 


4.2212147 (2) 

5.0080856 (2) 
5.9752941 (2) 
7. 1721622 (2) 
8.6638166 (2) 


1. 1431115 
1. 3974809 
1. 7197003 
2. 1310065 
2. 6603265 


(3) 
(3) 
(3) 
(3) 
(3) 


3. 1324761 
3. 9459207 
5. 0079025 
6. 4063588 
8. 2648601 


(3) 
(3) 
(3) 
(3) 

(3) 


8.6601663 (3) 
1.1240327 (4) 
1.4712160 (4) 
1.9428762 (4) 
2.5902033 (4) 


.71 
.72 
.73 
.74 

.75 


3.3716364 (2) 
3.9873716 (2) 
4.7463268 (2) 
5.6893286 (2) 
6.8711446 (2) 


1.0537046 (3) 
1.2908644 (3) 
1.5937432 (3) 
1.9841824 (3) 
2.4925823 (3) 


3. 3474663 
4. 2478017 
5. 439209K 
7. 03293S3 
9. 1892247 


(3) 
(3) 
(3) 
(3) 
(3) 


1. 0759010 
1. 4141363 
1. 8779832 
2. 5217787 
3. 4269672 


(4) 
(4) 
(4) 
(4) 
(4) 


3.4883257 (4) 
4.7489570 (4) 
6.5405516 (4) 
9.1209641 (4) 
1.2891309 (5) 


.76 

.77 
.78 
.79 
.80 


8.3660947 (2) 
1.0276297 (3) 
1.2743971 (3) 
1.5970122 (3) 
2.0243431 (3) 


3. 1618092 (3) 
4.0531621 (3) 
5.2556423 (3) 
6.9006308 (3) 
9. 1856068 (3) 


1. 2143196 
1. 6244726 
2. 2023403 
3. 0296111 
4. 2347605 


(4) 
(4) 
(4) 
(4) 
(4) 


4. 7175815 
6. 5857*93 
9. 3349272 
1. 3453030 
1. 9746667 


(4) 
(4) 
(4) 
(5) 
(5) 


1.8486494 (5) 
2.6930537 (5) 
3.9909042 (5) 
6.0258827 (5) 
9.2871234 (5) 


.81 
.82 
.83 
.84 
.85 


2.5985857 (3) 
3.3826250 (3) 
4.4722115 (3) 
6.0166633 (3) 
8.2551518 (3) 


1.2413348 (4) 
1.7058333 (4) 
2.3882318 (4) 
3.4141676 (4) 
4.9972061 (4) 


6. 0244685 
8. 7393424 
1. 2956045 
1. 9680620 
3. 0728282 


(4) 
(4) 
(5) 
(5) 
(5) 


2. 9572117 
4. 5284030 
7. 1085875 
1. 1473581 
1. 9109387 


(5) 
(5) 
(5) 
(6) 
(6) 


1.4640725 (6) 
2.3665859 (6) 
3.9336838 (6) 
6.7461970 (6) 
1.1985330 (7) 


.86 
.87 
.88 
.89 
.90 


1. 1582518 (4) 
1.6673365 (4) 
2.4726623 (4) 
3.7972735 (4) 
6.0789976 (4) 


7.5129608 (4) 
1. 1648160 (5) 
1.8715471 (5) 
3.1357024 (5) 
5.5223607 (5) 


4. 9500854 
8. 2655KW 
1. 4388026 
2. 6299573 
5. 0951320 


(5) 
(5) 
(6) 
(6) 
(6) 


3. 2984010 
5. 9314923 
1. 1186011 
2. 2300302 
4. 7538466 


(6) 
(6) 
(7) 
(7) 
(7) 


2.2165832 (7) 
4.2928225 (7) 
8.7706067 (7) 
1.9080277 (8) 
4.4730772 (8) 


.91 
.92 
.93 
.94 
.95 


1.0234053 (5) 
1.8334510 (5) 
3.5540054 (5) 
7.6378369 (5) 
1.8898335 (6) 


1.0330799 (6) 
2.0822980 (6) 
4.6133600 (6) 
1.1567756 (7) 
3.4349093 (7) 


1. 0591208 
2. 4017606 
6. 0816052 
1. 7791735 
6. 3399610 


(7) 
(7) 
(7) 
(8) 
(8) 


1. 0980174 

2. 8013155 
8. 1069617 
2. 7670704 
1. 1832720 


(8) 
(8) 
(8) 
(9) 
(10) 


1.1479956 (9) 
3.2950159 (9) 
1.0898261 (10) 
4.3398773 (10) 
2.2270714 (11) 


.96 
.97 
.98 
.99 
1.00 


5.7351988 (6) 
2.4035009 (7) 
1.8152204 (8) 
5.7774227 (9) 

CO 


1.3031098 (8) 
7.2819128 (8) 
8.2499417 (9) 
5.2518697 (11) 

CO 


3. 0066486 
2. 2402975 
3. 8073399 
4. 8476752 

CO 


(9) 

(10) 

(11) 

(13) 


7. 0146360 
6. 9691088 
1. 7700.103 
4. 5243605 


(10) 
(ID 
(13) 
(15) 


1.6503480 (12) 
2.1862541 (13) 
8.3603249 (14) 
4.2581255 (17) 

CO 



*The figures in parentheses indicate the power of ten by which the adjacent entry is to be multi- 
plied; e.g., 1.0412490 (2) =104.12490. 

Table 2. Simple power series (eq 2.1-4) 
These m-values are a lower bound on the number of terms to sum for &/(£) to ten significant figures. 





7=0 


;=8 




£2=0.01 


&2=0.99 


£2 = 0.01 


£2=0.99 


Minimum w-value for which 


3 


823 


4 


4615 
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Table 3. Difference sum (eq 2.2-6) 

The contribution of the difference sum to the total il,(k) is here compared with the contribution from the closed -form approximation. The m-values 
at which difference-sum terms contribute less than a part in 10" may be compared with m-values in table 2 for the simple power series. 



j=0 



n,-(A) (exact) 

Closed— form approximation: 

00 

]r}<a)^ 

771 = 

Difference — sum remainder: 

00 

^2[cx m (j)-a*Jj)]k^ 

W = 

Minimum w-value for which 
[a m (j)-al(j))k^ 
O/Gfc) 



3. 141651562 (0) 



7.071421389 (-1) 



2.434509423 (0) 



£2=0.01 



5.723176823 (0) 



2.825999794 (0) 



2.897177030 (0) 



3.147939532 (0) 



3.148259078 (0) 



-3. 195460509 ^-4) 



4.524360528 (15) 



4.524360528 (15) 



-4.282887696 (-4) 



Table 4. Bessel function for n 


lulation for 


Gj(A0 (eqs 2.3-27, 28) 




j= 


-0 


7 = 8 




£2=0.01 


£2 = 0.99 


fc2=0.01 


£2 = 0.99 


n,-(fc)- Exact 


3. 141 652 


5.723 177 


3. 147 940 


4.524 361X10+15 


rJ|[»(j4j)] 


3. 141 612 


3. 337 002 


3. 147 270 


1.512 709X10+ 3 


tt/o D 2 (;+20] 










-©KimHK^)] 










+i^(2j+l)/c[/c2(j+|)] 


3. 141 652 


3. 856 095 


3. 147 940 


1.527 585X10+* 


(3W64) (j+l) (7+g) £* 


1.2X10-9 


0.113 953 


1.2X10-7 


11.556 450 



Table 5. Error function formulation for ftj(k) {eq 2.4-41) 





j= 


= 


;=8 




£2 = 0.01 


£2=0.99 


£2 = 0.01 


£2=0.99 


«,-/ (k) -Exact 


3. 141 652 

3. 131 384 
3. 144 486 
3. 140 552 


5.723 177 

1. 781 883 

3. 121 919 

-1.354 046 


3. 147 940 

2. 994 801 
3. 177 138 
3. 159 417 


4.524 361X10+15 


VW(2;+l)/c 2 (l-A-2)-M-rf\ 
1st 2 terms in 2.4-41 
All 3 terms in 2.4-41 


3.822 991X10+U 
3.992 220X10+15 
3.958 982X10+15 



Table 6. Numerical integration technique (eqs 2.5, 44~49) 





j= 


=0 


7 = 8 


N 


£2=0.01 


£2=0.99 


£2=0.01 


£2=0.99 








fly(£) 




2 
3 

4 
5 
6 


3. 141 632 
3. 141 644 
3. 141 653 
3. 141 652 
3. 141 652 


7. 701 552 
6. 826 700 
5.931 746 
5. 794 716 
5. 615 818 


3. 145 825 
3. 147 146 
3. 148 080 
3. 148 015 
3. 147 932 


52.359 878X10+15 
39.269 908X10+15 
24.434 610X10+'5 
20.725 792X10+15 
15. 333 964X10+15 


10 


3. 141 652 


5.601 518 


3. 147 940 


8.430 212X10+1* 


oo (Exact) 


3. 141 652 


5.723 177 


3. 147 940 


4. 524 361X10+15 
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Table 2, which applies to what has been called the simple power series above, illustrates 
how rapidly the number of terms required increases as either k or j goes up. These figures 
were obtained by setting the m'th term in eq. (2.1-4) equal to fy(&)X10 -11 and solving for m. 5 
For large j and k values, the use of this formulation becomes almost prohibitively laborious. 

In table 3, where the difference sum technique is examined, a number of interesting charac- 
teristics of this formulation are apparent. For small values of k, the improvement over the 
simple power series is negligibly small for all values of j. For k approaching unity, there is a 
considerable reduction in the labor involved to obtain any predetermined accuracy. For 
j=0, k 2 =0.99 to obtain a maximum of ten significant figures, note that the number of terms 
required is cut from 823 to 523. While this improvement is not as spectacular as some cases 
which have been examined [8], it is by no means unimportant. The real gain in using this 
technique is seen at i=8, P=0.99; here we see that the closed-form approximation alone gives 
tij(k) correct to one part in 10 19 . This can be compared with the more than 4615 terms indicated 
in table 2 for ten or less significant figures using the unmodified power series. 

To check the behavior of the Bessel function formulation, section 2.3 above, the BAAS 
tables [16] of I (x) and Ii(x) were used along with suitable interpolation tables [14]. The 
results are given in table 4. It will be recalled that these relations were derived essentially 
by approximating only the first few terms of the integrand in eq 2-1. For this reason, the 
formulation as would be expected yields excellent values of Slj(k) when P=0.01 at both ^'=0 
and j=S) but for k 2 near unity, this technique of calculation is entirely inadequate. 

The error function formulation [17], section 2.4, table 5, is also based on approximations 
to only the first few terms in the expansion of the integrand. Note that the leading term from 
2.4-41 is not a bad approximation to the value of Qj(k); but the inclusion of the second or 
third terms results in little or no improvement. Evidently many more terms are needed to 
obtain useful results by this technique. 

The importance of the erf and I n formulations of Qj(k) arises, then, not so much for their 
usefulness in numerical computation, as from the relations they demonstrate between the new 
transcendental function considered here and older, more familiar forms. 

The numerical integration formulas, from section 2.4, are quite different in kind. These, 
it will be observed from table 6, give excellent agreement for small k values, although they 
leave something to be desired for k near unity. It is interesting to note that, at these selected 
points, the accuracy of the formulas does not increase uniformly as N increases — for example, 
at k 2 =0.99, j=0, the formulas with iV=5 and 6 are virtually equivalent, and each is better 
than A 7 =10. The approximation for large jk 2 is excellent. Note that by substituting the 
value of Uj(k) at k 2 =0.99 and .7 = 8 from table 1 it may be concluded that 

lim w N (0)l^w N (i) = lim |=^J f* (<-l)(*-2)(*-3) ... (t-N)dt s 0.0144 . . . 

iV->co / j=0 A^co 1\ • iV! Jo 

It thus appears that there is a multiplicity of choices available for computing Qj(k) for 
small k values, and for k approaching unity the recursion formula (1-7) can be used with the 
appropriate high-& expansions (1-10, 11) for ordinary elliptic integrals. However, for high- 
precision coverage of the entire range 0<&<O with a single formulation, only the difference- 
sum technique appears to be practical. 



5 It should be noted that, since (2.1-4) is not an alternating series, the solution for m of the inequality 

amO"H- 4m <fl/(fc)X10-" 
serves only as a general indication of the convergence rate and does not necessarily imply that the remainder 

00 

ra + i 

is less than fi,-(fc) X10 -11 . It does, however, indicate the minimum number of terms to sum, beyond which there exists, for A: 2 <1, a value m'>rn for 
which 

Rm f> m(j, 
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The authors thank K. J. Herbold for programing and computing table 1, R. L. Bach for 
numerical checks by alternative formulations, L. C. Maximon for suggesting Stirling's formula 
as a basis for the difference-sum, and M. J. Berger for other helpful suggestions in the comple- 
tion of (1) is work. 

4. Appendix A. Evaluation of a* 0") 
Stirling's formula [11, p. 5] may be written in the form 

z\=^2^z z e- z H{z) (Al) 

in which H(z) can be represented by the asymptotic expansion [18] 

limi?(^)-l+l/(12^) + l/(288^ 2 )-139/(518402 3 )- .... (A2) 

In order to use formula (Al) above, it is convenient to write eq (2-5) as 

/ :\ T (4m)! J. / 4m , A /AoN 

a '" U) = (64T (2«)l(ml)' Ji 0^T +1 ) (A3) 

in which the product is to be taken as unity for ,7 = 0. Substitution of (Al) for the factorials 
in (A3) gives 

*-CA=(l/mV2) H Z)mm) i (^l + O- (A4) 

By using (A2J for the //terms, and expressing the finite product as a /th degree polynomial 
in m, it follows that 






A 1 



f2x— 1 



where 



^ (2j)! 96 JU ; { J J j+ (2j)! 4^ + (2j)lJ 

= (lA2)EC0')m'-'-' (A5) 



f'o(i) = +2+ 3 ^ (A6-1) 

C(i) = +2 +3 ^-^(4i 2 -3) (A6-2) 

a(j) = +2+ 3 ^ — (48j^-64j :, -72.r+16i+27) (AC-3) 



and also 



<%")=2- 3 '( 2 £) (A6-4) 
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which can be established by methods of mathematical induction. Now it is possible to show 
that m n can be expressed in terms of binomial coefficients, in the form 



m 



k=o\ m / 



where the jft^s are related to Kummer's numbers [19]. It can readily be shown that 

0n,o=(— l) n and p n ,n=n\ 
A few values of these constants are tabulated below 



V 

\ 





1 


2 


3 


F 

4 



1 


+ 1 


-1 

+ 1 


+ 1 
-3 

+ 2 


-1 

+ 7 
-12 

+ 6 


+ 1 
— 15 


2 




+ 50 


3 






— 60 


4 








+ 24 















(A7) 



Now ignoring terms of the order of 1/m 2 in the expression (A-5), for large m values, a m (j) 
can be approximated by the expression 



(1/V2) [g tf^m'-'-'+^COm- 1 ]. 



Replacing 1/m by its equivalent to this order, l/(m+l), and using (A-7) in the sum term, it 
follows that 

«£(i)=(i/V2) ["S C t (j) 's &_,_,.* ( m + k )+C,(j) -q-r"| 

U=o k=o \ m / m-f-ij 

■w j ' CO {^+5 B - ( m+ t <_1 )} <A -« 

which is eq 2.2-8 of the text above. The B tj expression in eq (A-8) is obtained by combining 
the C t (j) and p ntk values as indicated and factoring out Cj(j). The expression (2.2-11) above 
is then obtained using the familiar series 



and 



In (l-s) = — Z) x m+1 /(m + l) (A9) 

m=0 

(i _*)-*-!= ±( m + k )ir. (Aio) 
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